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Abstract 

The direct CP violation parameter Re(e'/e) has been measured from the decay 
rates of neutral kaons into two pions using the NA48 detector at the CERN SPS. 
The 2001 running period was devoted to collecting additional data under varied 
conditions compared to earlier years (1997-99). The new data yield the result: 
Re(e'/e)=(13.7 ± 3.1) x 10~ 4 . Combining this result with that published from the 
1997, 98 and 99 data, an overall value of Re(e'/e)=(14.7 ± 2.2) x 10~ 4 is obtained 
from the NA48 experiment. 

(To be published in Physics Letters B) 
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1 Introduction 

CP violation was discovered 38 years ago in the decays of neutral kaons flU . Recently 
CP violation in the B° mesons has also been observed 0; nevertheless neutral kaons 
remain a privileged system for the study of the phenomenon. 

CP conservation would imply that the Kg and Kl particles are pure CP-eigenstates 
and that Kl decay only into CP=— 1 and K$ into CP=+1 final states. The observed 
signal of the forbidden Kl — > 27r decays (CP=+1) indicates that CP is not a conserved 
symmetry. 

CP violation can occur via the mixing of CP eigenstates, called indirect CP viola- 
tion, represented by the parameter e. CP violation can also occur in the decay process 
through the interference of amplitudes with different isospins. This is represented by the 
parameter e' and is called direct CP violation. 

In the Standard Model of electro-weak interaction, CP violation is naturally ac- 
commodated by an irreducible complex phase in the quark mixing- matrix ||. Current 
theoretical predictions of e'/e range ~— lOxlCr 4 to ~+30xl0~ 4 (4|. 

Experimentally, it is convenient to measure the double ratio R of decay widths, 
which is related to the ratio e'/e as follows: 

T(K L -> 7r°7r°) T(K L -> tt+tt-) , 
R = -A— - —1 — i ss 1 - 6 x Re e'/e 1 

T(K S -> 7l°7l°) T(K S -> 7T + 7T-) V ' ' V ' 

In 1993, two experiments published their final results: NA31 measured Re(e'/e) = 
(23.0 ± 6.5) xl(T 4 , and the result of E731 § was Re(e'/e) = (7.4 ± 5.9)xl0~ 4 . Recently, 
two experiments announced results from samples of their total statistics. NA48 published 
a result of Re(e'/e) = (15.3 ± 2.6) xl(T 4 , using data collected in 1997 0, 1998 and 99 §, 
and KTeV presented a preliminary result of Re(e'/e) = (20.7 ± 2.8)xl0~ 4 || on data 
accumulated in 1996 |L0| and 97. These observations confirmed the existence of a direct 
CP- violation component. 

This paper reports a measurement of Re(e'/ e) performed using the 2001 data sample, 
recorded in somewhat different experimental conditions by the NA48 experiment. 

After the 1999 data-taking period, the drift chambers of the experiment were dam- 
aged by the implosion of the beam tube. The data taking in 2001 took place with rebuilt 
drift chambers. Thanks to the possibility of a better SPS duty cycle, increased by a factor 
1.8 with respect to the 1998-99 running period, the data could be taken at a 30% lower 
beam intensity, allowing the insensitivity of the result to intensity-related effects to be 
checked, and the statistics for the final e'/e measurement by NA48 to be completed. The 
statistics accumulated during the 93 days of the 2001 data-taking period is roughly half 
of the total statistics accumulated in the 263 days of the 1998 and 99 periods. 

Details of the apparatus and of the data analysis can be found in || , here only the 
differences with respect to the 1998-99 data-taking will be stressed. 



2 The method 

Re(e'/e) is derived from the double ratio R. The experiment is designed to exploit 
cancellations of systematic effects contributing symmetrically to different components of 
the double ratio. 

The four decay modes are collected simultaneously, which minimises the sensitivity 
of the measurement to accidental activity and to variations in beam intensity and detection 
efficiency In the analysis K$ events are further weighted by the Kl/Ks intensity ratio to 
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eliminate the small and slow variations of the Kl and K$ beam intensities. To maintain 
the simultaneous data-taking of 7r°7r° and 7r + 7r _ decays, dead-time conditions affecting 
one mode are recorded and applied offline in all modes. 

Kl and K$ decays are provided by two nearly-collinear beams with similar mo- 
mentum spectra, converging to the centre of the main detector. The same decay region is 
used for all modes. In order to minimise the acceptance correction due to the difference in 
mean decay lengths, Kl decays are weighted as a function of their proper lifetime, such 
that the Kl decay distribution becomes similar to that of K$- In this way, the accuracy of 
the result does not rely on a detailed Monte Carlo simulation of the experiment and only 
small remaining differences in beam divergences and geometries need to be corrected using 
Monte Carlo simulation. To be insensitive to residual differences in the beam momentum 
spectra, the analysis is performed in bins of kaon energy. 

Ks decays are distinguished from Kl decays by a coincidence between the decay 
time and the registered times of the protons producing the Ks beam. As the same method 
is used for 7r + n~ and 7r°7r° decays, the double ratio is sensitive only to differences in 
misidentification probabilities between the two decay modes and not to their absolute 
values. 

Finally, high-resolution detectors are used to detect the 7r + 7r~ and 7r°7r° final states 
in order to minimise residual backgrounds which do not cancel in the double ratio. 



3 Beams and detectors 
3.1 Beams 

The Kl and Ks beams [1 1] are produced in two different targets by protons from 
the same CERN SPS beam. In the 2001 run the SPS had a cycle time of 16.8 s with a spill 
length of 5.2 s and a proton momentum of 400 GeV/c q Since the Ks and Kl beams are 
produced concurrently, the Ks/Kl ratio is maintained stable to within ±10%. 

The primary, high-flux proton beam (~2.4xl0 12 protons per pulse) impinges on 
the Kl target (a 400 mm long, 2 mm diameter rod of beryllium), with an incidence 
angle of 2.4 mrad relative to the Kl beam axis. The charged component of the outgoing 
particles is swept away by bending magnets. The neutral beam passes through three 
stages of collimation and the fiducial region starts at the exit of the "final" collimator, 
126 m downstream of the target. At this point, this neutral beam is dominated by long- 
lived kaons, neutrons and photons; only small fractions of the most energetic short-lived 
components (Ks and A) survive. 

The protons not interacting in the Kl target are directed onto a mechanically 



bent mono-crystal of silicon [|12|| . A small fraction of protons satisfies the conditions for 
channelling and is bent to produce a collimated beam of ~5xl0 7 protons per pulse, 
which is then deflected back onto the Kl beam axis and finally directed to the Ks target 
(of similar dimensions as the Kl target) located 72 mm above the Kl beam axis. A 
combination of a sweeping magnet and collimator selects a neutral beam at 4.2 mrad to 
the incoming protons. The decay spectrum of kaons at the exit of the collimator (6 m 
downstream of the target) is similar to that in the Kl beam, with an average energy of 
110 GeVQ. Two -pion decays from this beam come almost exclusively from Ks decays. 



^ The cycle time was 14.4 s, the spill length 2.38 s and the proton momentum 450 GeV/c in the 1998 
and 99 runs. The effective spill length, given by the remaining time structures in the beam is «3.6 s 
in 2001, compared to sal. 7 s for the 1998-99 data. 
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The tagging station (Tagger) is located on the path of the K$ proton beam after 
the bent crystal. It consists of two ladders of 12 scintillator strips each, covering the beam 
horizontally and vertically pT3 1 . A proton crosses at least two scintillators, one horizontal 



and one vertical. The reconstructed time per counter has a resolution of ~140 ps, and 
two pulses 4-5 ns apart can be resolved. 

The beginning of the Ks decay region is sharply defined by an anti-counter (AKS), 



located at the exit of the K$ collimator ||14|| . It is composed of a photon converter followed 
by three scintillator counters. Its main purpose is to veto all upstream decays from the 
Ks beam. 

The decay region is contained in an evacuated (< 3 x 1CT 5 mbar) 90 m long tank 
with a 0.9 mm (0.003 radiation length) thick polyamide (Kevlar) composite window at 
the end. The neutral beam continues in a 16 cm diameter evacuated tube to the beam 
dump, downstream of all detector elements. 



com- 



3.2 Detectors 

Charged particles from decays are measured by a magnetic spectrometer |15 
posed of four drift chambers with a dipole magnet (inducing a transverse momentum-kick 
of 265 MeV/c in the horizontal plane) between the second and third chambers. These 
chambers and their interconnecting beam tube are aligned along the bisector between the 
converging Kg and Kl beam axes. Each chamber is comprised of eight planes of sense 
wires, two horizontal, two vertical and two along each of the 45° directions. In the third 
chamber, only the horizontal and vertical planes are instrumented. The average efficiency 
per plane is 99.5%, with a radial uniformity better than ±0.2%. The space point reso- 
lution is ~ 95 /im. The momentum resolution is cr(p)/p = 0.48% © 0.009% x p, where 
the momentum p is in GeV/c. These performance figures are similar to those obtained 
previously |§. 

The magnetic spectrometer is followed by a scintillator hodoscope, composed of two 
planes segmented in horizontal and vertical strips. Fast logic combines the strip signals 
(arranged in four quadrants) for use in the first level of the ir + n~ trigger. 

A liquid Krypton calorimeter (LKr) is used to reconstruct 2tt° decays. It is 
a quasi- homogeneous detector with an active volume of ~10 m 3 of liquid Krypton. Cu- 
Be-Co ribbon electrodes of 40 /j,m x 18 mm x 125 cm define ~ 13000 cells (each with 
2 cm x 2 cm cross-section) in a structure of longitudinal projective towers pointing to 



the centre of the decay region [16]. The calorimeter is 27 radiation lengths long and fully 
contains electro- magnetic showers with energies up to 100 GeV. The energy resolution of 
the calorimeter is cr(E)/E = (3.2 ± 0.2)%/ © (9 ± 1)%/E © (0.42 ± 0.05)% with E in 



GeV [|l7fl. The LKr calorimeter also is used, together with an iron-scintillator calorimeter, 
to measure the total deposited energy for triggering purposes. 

Finally, at the end of the beam line, muon counters are used to identify K L —> -k\w 
(K /t3 ) decays. 

Two beam counters are used to measure the intensity of the beams: one is located at 
the extreme end of the Kl beam line {Kl monitor) and the other (Kg monitor) views the 
Kg target station. For the 2001 data taking, another Kl monitor with a higher counting 
rate was added and a Ks monitor near the tagging station was installed. These allow 
better measurements of the beam structures to be made down to a time scale of ~ 200 ns. 



2 ) Despite the different proton momentum in the 1998 and 99 runs, the Kl and Ks spectra remain 
similar, depending only on the choice of production angles to compensate for the length of the Ks 
collimator. 
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3.3 Triggers 

The rate of particles reaching the detector is around 400 kHz. The trigger is designed 
to reduce this rate to less than 10 kHz, with minimal loss from dead time and inefficiencies. 
A part of the read-out rate is reserved for redundant low-bias triggers that collect data 
used for the direct determination of the trigger inefficiencies. 

Triggers initiated by the beam monitors are used to record the accidental activity, 
with rates proportional to and K$ decay rates. Beam monitor signals are down-scaled 
and delayed by 69 fis corresponding to the periodicity of the slow proton extraction (3 
SPS revolutions). 

3.3.1 Trigger for 7i°7i° decays 

The trigger for tt°tt decays [IS| operates digitally on the analogue sums of signals 



from 2x8 cells (in both horizontal and vertical orientations) of the LKr calorimeter. 
These sums are converted into kinematic quantities by a "look-up table" system. 

The trigger requires an electro-magnetic energy deposit greater than 50 GeV, a 
centre of energy (distance between the extrapolated kaon impact point at the calorimeter 
plane and the beam axis) smaller than 25 cm, and a decay vertex less than 5 K$ lifetimes 
(ts) from the beginning of the decay volume. Requesting less than 6 peaks within 9 ns in 
both projections helps to reject background from Kl — > 37r°(this condition is released if 
accidental activity is detected close in time). 

A trigger for 3n° decays, given by the down-scaled 7r°7r° trigger without the peak 
condition, is used for tagging studies. 



3.3.2 Trigger for 7T + 7T~ decays 

The 7r + 7r~ decays are triggered with a two-level trigger system. At the first level, 
the rate is reduced to 100 kHz by a coincidence of three fast signals: opposite quadrant 
coincidence in the scintillator hodoscope (Q x ), hit multiplicity in the first drift chamber 
integrated over 200 ns, and the total calorimetric energy (E tot , with a threshold of 35 GeV). 

The second level of the 7r + 7r~ trigger consisting of hardware coordinate builders 
and a farm of asynchronous microprocessors, reconstructs tracks using data from the drift 
chambers. Triggers are accepted if the tracks converge to within 5 cm, their opening angle 
is smaller than 15 mrad, the reconstructed proper decay time is smaller than 4.5 t$ and 
the reconstructed inr mass is greater than 0.95 my. 



4 Event reconstruction and selection 
4.1 7r°7r° 

K — > 7r°7T decays are selected using only data from the LKr calorimeter. The re- 
construction of photon showers and the details of the small corrections applied to the 
energy and position measurements can be found in || and ||20|| . Photon showers in the 
energy range 3 — 100 GeV are used. Fiducial cuts are applied to ensure that the photon 
energies are well measured: the shower position should be more than 15 cm away from 
the axis of the beam tube, more than 11 cm away from the outer edges of the calorimeter 
and more than 2 cm away from a defective calorimeter channel (~ 0.4% of the channels). 
7r°7T decays are selected by requiring four showers which are reconstructed within ± 5 ns 
of their average time and fulfill the cuts above. The minimum distance between photons 
is required to be more than 10 cm. To reduce the background from Kl — > 3tt° decays, 
events are rejected which have an additional cluster of energy above 1.5 GeV and within 
± 3 ns of the time of the 7r°7r° candidate. 
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From the measured photon energies Ei and impact point positions on the calorime- 
ter Xi,i/i, the distance D from the decay vertex to the calorimeter is computed as follows, 
assuming that the invariant mass of the four showers is the kaon mass (m^): 



D = JY li Y lj>i E i E j ({x i - Xj) 2 + (yi - yj) 2 )/m K - 



(2) 



The average resolution on the decay vertex position is about 55 cm, and the resolution 
on the kaon energy is ~ 0.5%. 

The invariant masses m 1 and m 2 of the two photon pairs are computed using D and 
compared to the nominal 7r° mass (m^o). For this, the following x 2 variable is constructed: 



X 



(mi + m 2 )/2 - m^o 



+ 



(mi - m 2 )/2 



(7_ 



(3) 



The mass combinations mi + m 2 and m x — m 2 are to good approximation uncorrelated. 
a + and er_ are the resolutions of (mi + m 2 )/2 and (mi — m 2 )/2 parameterised from the 
data as a function of the lowest photon energy. Typical values of a + and <r_ are 0.4 and 
0.8 MeV/c 2 . Out of the three possible photon pairings, the one with the lowest x 2 value is 
kept. To select good 7r°7r° candidates and reject the residual background from K L — > 3ir° 
decays, the cut x 2 < 13.5 is applied. 

The event time is computed by combining eight time estimators from the two most 
energetic cells of each cluster. An average resolution of 220 ps is thereby obtained. 



4.2 7T+7T- 

The n + n~ events are reconstructed from tracks using hits in the drift chambers of 
the spectrometer; the track momenta are calculated using the measured magnetic field 
map and alignment constants. 

A vertex position is calculated for each pair of tracks with opposite charge after 
correcting for the small residual magnetic field due to the magnetisation of the vacuum 
tank (~ 2xl0~ 3 Tm). The average resolution on the longitudinal vertex position is about 
50 cm, whereas the transverse position resolution is around 2 mm. Since the beams are 
separated vertically by about 6 cm in the decay region, a clean identification of K$ and 
Kl decays is possible using the reconstructed vertex position. 

Only tracks with momenta greater than 10 GeV/c and not closer than 12 cm to 
the centre of each chamber are accepted. The separation of the two tracks at their closest 
approach is required to be less than 3 cm. The track positions, extrapolated downstream, 
are required to be within the acceptance of the LKr calorimeter and of the muon veto 
system, in order to ensure proper electron and muon identification. 

The kaon energy is computed from the opening angle 9 of the two tracks upstream 
of the magnet and from the ratio of their momenta p\ and p 2 , assuming a K — > 7r + 7r~ 
decay: 

E K = M(m 2 K -nml) where ft = ^ + ^ + 2 (4) 
V 9 2 p 2 pi 

This measurement of the kaon energy is independent of the absolute magnetic field and 
relies mostly on the knowledge of the geometry of the detector. 

A variable A related to the decay orientation in the kaon rest frame is defined as 
A = |pi -p 2 \/{pi +P2)- A cut is applied to A {A < mm(0.62, 1.08-0.0052 x E K ), where 
E K is in GeV), to remove asymmetric decays in which one of the tracks could be close to 



5 



the beam tube where the Monte Carlo modelling is more critical. This cut also removes 
A — ► pn~ decays. 

To reject background from semileptonic Kl decays, events with tracks consistent 
with being either an electron or a muon are rejected. To identify electrons, the ratio E/p 
of the energy of the matching cluster in the LKr calorimeter to the track momentum is 
computed. Pion candidate tracks are required to satisfy E/p < 0.8. Tracks are identified 
as muons if hits are found in time in the muon counters near the extrapolated track impact 
point. 

For good 7T + 7r~ events, the reconstructed mass m nn should be consistent with the 
kaon mass. The resolution on the invariant mass a m is typically 2.5 MeV/c 2 . An energy- 
dependent cut at ±3<T m is applied. A further reduction of background from semileptonic 
decays is achieved with a cut based on the transverse momentum of the kaon. To define a 
selection which is as symmetric as possible between Kl and Ks decays, the variable pt is 
used, defined as the component of the kaon momentum orthogonal to the line joining the 
production target and the point where the kaon trajectory crosses the plane of the first 
drift chamber. To select 7r + 7r~ candidates, the cut pt' 2 < 2 x 10" 4 GeV 2 /c 2 is applied. 

The time of the n + ir~ decay is determined from hits in the scintillator hodoscope 
associated with the tracks and has a resolution of ~150 ps. The events with insufficient 
information to determine the decay time accurately are discarded. This inefficiency is 
0.1% and is measured to be equal for Ks and K L . 



4.3 K s tagging 

A decay is labelled Ks if a coincidence is found (within a ±2 ns interval) between 
its event time and a proton time measured by the Tagger. Figure [I] shows the time 
distributions for Ks and Kl decays to 7r + 7r~ which have been identified as such by their 
vertex positions in the vertical direction. A similar procedure is not possible for 7r°7r° 
decays; therefore tagging provides the only way to distinguish Ks from Kl- The selection 
of Ks and Kl samples is done by means of tagging for both the 7r + 7r~ and 7r°7r° modes, 



K L — >7l ' K 




Time(K )-time(proton) (ns) 



Figure 1: Time coincidence for Ks and Kl ^ + tt decays, identified by their reconstructed 
vertex. 
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so that systematic effects are mostly symmetric. 

The probability that a Kg decay is assigned to the Kl beam, due to coincidence 
inefficiencies, is denoted by a$L- This tagging inefficiency is directly measured in the 
7T + 7r~ mode, using the vertical vertex position, and is (1.12 ± 0.03) xlO -4 , dominated by 
the Tagger inefficiency. The inefficiency could be different for 7r + 7r~ and decays, since 
different detectors are used. The difference between the two modes is estimated using a 
large sample of Ks and Kl decays into 27r° and 3tt° where one of the photons converts 
into an electron-positron pair. The time from the LKr clusters can then be compared with 
the scintillator-hodoscope time from the two tracks. The conclusion is that the 7r°7r° and 
7T + 7r~ tagging inefficiencies agree within an uncertainty of ±0.5xl0 -4 , which corresponds 
to an uncertainty on R of ±3xl0 -4 . The tagging inefficiency can also be measured directly 
in 7T 7r° events with a subsequent 7r° — > e + e~7 Dalitz decay, which allows the Ks to be 
identified from the two-track vertex. However this method is statistically limited. Another 
cross-check is performed in special runs where only Ks are present; it gives a tagging 
inefficiency in agreement with the result above. 

The probability that a Kl decay is identified as a Ks decay due to an accidental 
coincidence between the event and a proton is called accidental tagging and is denoted as 
otLs- It is measured in the 7r + 7r~ mode to be (8.115 ± 0.010) x 10™ 2 Q. The c^ls difference, 
Ao>ls, between the 7r 7r° and 7r + 7r~ modes is estimated by measuring the probability to find 
a proton within time windows 4 ns wide, located before or after the event time in tagged 
Kl events (i.e. events with no proton in coincidence). Ten windows are chosen, centred 
at 5 ns intervals from the event time so as to follow the 200 MHz structure of the proton 
beam. The extrapolation from the side windows to the coincidence window is performed 
in the 7t + tc~ mode with vertex-selected Kl and in the 7r°7r° mode using 3n° decays, which 
come almost entirely from the Kl beam (the very small contribution from the Ks beam 
being subtracted). The measured value of Aa^s = a^s 00 — o.ls ± = (+3.4 ± 1.4) xl0~ 4 
corresponds to a correction on R of (+6.9±2.8) x 10~ 4 . The origin of this effect is discussed 
in section |6.2| . 

4.4 Definition of the decay region 

The fiducial ranges in kaon energy Ek and in proper time r used to count events 
are chosen to be 70 < E K < 170 GeV and < r < 3.5 r s , where r = is defined at the 
position of the AKS counter and Ts is the Ks mean lifetime. For Kl events, the decay 
time cut is applied on reconstructed r, while for Ks events the cut at r = is applied 
using the AKS to veto decays occurring upstream. The nominal r = positions defined by 
the AKS differ by 21.0 ± 0.5 mm between 7r + 7r~ and 7r°7r° decays. The veto inefficiency is 
0.36% for 7r°7r° events and 0.22% for 7r + 7r~ decays. Given the fractions of decays occurring 
upstream of the AKS position (respectively 5.8% and 4.0%), the correction to be applied 
to the double ratio is (1.2 ± 0.3) x 10" 4 . 

Ks — > 7T7T decays can be produced in both beams by scattering of beam particles 
in the collimators and, in the case of the Ks beam, in the AKS counter. To reduce 
this contamination, a cut on the extrapolated kaon impact point at the level of the LKr 
calorimeter is applied to all events. For 7r°7r° decays, it is defined as the energy-weighted 
average x, y positions of the four showers at the face of the LKr calorimeter. For 7r + 7r~ 
decays, it is the momentum-weighted average position of the tracks measured upstream 

3 > It was (10.649 ± 0.008) x 10~ 2 in the 1998-99 data sample. 
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of the spectrometer magnet and projected onto the face of the LKr calorimeter. The 
extrapolated kaon impact point is required to be within 10 cm of the intersecting beam 
axes. The radii of the Kl and Kg beam spots are respectively 3.6 cm and 4.6 cm, so 
effects related to resolution smearing are negligible. 

4.5 Data quality selection 

Data collected in the mode are affected by an overflow condition in the drift 

chambers which resets the front-end readout buffers when there are more than seven hits 
in a plane within a 100 ns time interval. This occurs mostly when an accidental particle 
generates an electromagnetic shower upstream of the spectrometer and sprays the drift 
chambers with particles. To maintain the highest reconstruction and trigger efficiencies, 
events used in the analysis are required to have no overflows within ± 312 ns of the event 
time. To minimise possible effects of Kl/Ks intensity variation and to equalise the beam 
intensities seen by good events, this requirement is applied to both 7r + 7r~ and decays. 
The resulting event loss is 11% in the 2001 data sample. This is significantly smaller than 
the 20% loss observed in the 1998-99 data, due to the lower instantaneous beam intensity 
and to smaller noise in the drift chambers. Other dead time conditions affecting the first 
and the second level 7r + 7r~ trigger (« 0.3% 0) are also recorded and applied in the analysis 
to 7T 7r° decays. 

Because of large beam intensity variations at the very beginning of the spill, data 
from the first 0.2 s are not used in the analysis. The corresponding loss of events is ~ 1% 
and cancels in the ratio between 7r°7r° and 7r + 7r~ decays. 

5 R corrections and systematic uncertainties 

The data are divided into 20 bins in kaon energy, each 5 GeV wide. The numbers of 
K$ and Kl candidates are corrected for the mistagging probabilities discussed previously. 
The total numbers of events are 1.546 xlO 6 K L -> 2.159x 10 6 K s -> 7r°7r°, 7.136 x 10 6 

K L -> 7T+7T- and 9.605 x 10 6 K s -> tt+tt-. 

Corrections for trigger efficiencies, background subtractions and residual acceptance 
differences between Kl and Ks are applied separately in each energy bin before computing 
the average of R. 

5.1 Trigger efficiencies 

The 7r°7r° trigger efficiency is measured using a control sample of events triggered 
by a scintillating fibre detector located inside the LKr calorimeter. The efficiency is found 
to be (99.901 ± 0.015)%. The small inefficiency is Ks-Kl symmetric and no correction to 
the double ratio need be applied. 

The 7r + 7r~ trigger efficiency is (98.697 ± 0.017)% Q. The difference between the 
7T + 7r~ trigger efficiency for Ks and Kl decays is computed in each energy bin. The overall 
correction on the double ratio is (5.2 ± 3.6) x 10 -4 , where the uncertainty is given by the 
statistics of the control samples used to measure the efficiency. 

5.2 Backgrounds 

4 > It was 1.6% in the 1998-99 data. 

5 - 1 It was (97.782 ± 0.021)% in 1998-99. The improvement comes from the lower beam intensity and the 
better efficiency of the drift chambers in 2001. 
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5.2.1 Background to the 7r°7r° mode 

The background to the Kl — > tt 7i° signal comes uniquely from Kl — ► 37r° decays, 
while the Ks mode is background free. The Kl — > 37r° background has a flat x 2 distri- 
bution. To estimate this background, a control region is defined by 36 < x 2 < 135. The 
excess of Kl candidates in this region over a Monte Carlo expectation for -k°7t° decays is 
used to extrapolate the background in the signal region. The Monte Carlo includes the 
effect of non-Gaussian tails in the calorimeter resolution so as to reproduce the observed 
distribution in the Ks sample. 

The background is subtracted from the Kl — > 2tt° sample in bins of kaon energy 
and the resulting correction on the double ratio, taking into account the uncertainties in 
the non-Gaussian tails and in the background extrapolation, is (—5.6 ± 2.0) xlO -4 . 

5.2.2 Background to the it + it~ mode 

The background from A — > pir~ in the Ks — > 7r + 7r~ sample is negligible after the 
cut on the track momentum asymmetry. 

The residual K e 3 and K m3 backgrounds in the Kl sample are estimated by defining 
two control regions in the m n7T -pT' 2 plane. The first region, 9.5 < (m^ — m^) < 19.0 
MeV/c 2 and 300 < p T ' 2 < 2000 MeV 2 /c 2 , is dominated by K e3 events, while the second, 
-17.0 < (m^ - m K ) < -12.0 MeV/c 2 and 300 < p T ' 2 < 500 MeV 2 /c 2 , contains roughly 
equal numbers of K e3 and K m3 events. 

The background distributions in the control regions are modelled by a K e3 sample, 
selected with E/p > 0.95, and by a K m3 sample, obtained by reversing the muon veto 
requirement; the tails in the Kl — > 7r + 7r~ distribution are estimated from the Ks sample. 
The result is then extrapolated to the signal region. 

The overall K e3 background fraction is 10.5 xlO -4 , the K^ 3 background is 4.0 xlO -4 . 
The background subtraction is applied in bins of kaon energy and the resulting correction 
on the double ratio is (14.2 ± 3.0) xlO -4 , where the error has been estimated by changing 
the control regions and the modelling of the resolution tails. 

Kaon decays to 7r + 7r~7 have been shown to have a negligible effect on R [|J. 

5.2.3 Collimator scattering 

In the Ks beam, the cut on the extrapolated kaon impact point is stronger than 
the pt' 2 cut applied to 7t + tt~ decays, and therefore the contribution of beam scattering is 
removed symmetrically from both final states. On the contrary, in the Kl beam, the px' 2 
cut which is applied only in the 7r + 7r~ mode is stronger and therefore the small residual 
contribution from scattered events must be subtracted from the 7r°7r° sample. 

The correction for this asymmetry is computed from reconstructed Kl — * 7r + 7r~ 
candidates with an inverted p^ 2 cut. The scattered events are extracted from the peak 
at the kaon mass in the m ww invariant mass distribution. The correction to R is applied 
in bins of energy and it amounts to (—8.8 ± 2.0) xl0~ 4 . 

5.3 Acceptance 

The Ks and Kl acceptances are made very similar in both modes by weighting 
Kl events according to their proper decay time. The weighting factor takes into account 
the small interference term. A small difference in acceptances remains, related to the 
differences in Ks and Kl beam sizes and directions. This residual correction is computed 
using a large-statistics Monte Carlo simulation (4xl0 8 generated kaon decays per mode). 
The largest contribution to the correction comes from the difference between the Ks and 
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Kl beams near the beam axes in the spectrometer for 7r + 7r~ decays. The acceptance 
correction related to the tt°tt° mode is small. 

The systematic uncertainty on the acceptance correction is evaluated by varying 
the Ks beam halo, the beam positions and shapes, and the drift-chamber inefficiencies. 
The resulting systematic uncertainty is ±3.0 x 10~ 4 . A detailed comparison between a fast 
simulation and a GEANT [^TJ based simulation of the spectrometer was performed on the 
1998-99 sample. This resulted in an additional systematic uncertainty of ±2.3 x 10~ 4 . The 
final correction to R for the acceptance is: A-R(acceptance) = (+21.9 ± 3.5(MCstat) ± 
4.0(syst))xl0~ 4 . 



5.4 Energy and distance scales 

The determinations of the kaon energy, the decay vertex and the proper time in the 
7T 7r° mode rely on measurements of the photon energies and positions with the calorime- 
ter. 

The absolute energy scale is adjusted using Ks — > 7r 7r° decays. The energy scale 
is set such that the average value of the reconstructed decay position in a range centred 
around the anticounter matches the value found in a Monte Carlo simulation. This mea- 
surement of the energy scale is checked using data taken during special runs (so-called 
rj runs) with a tt~ beam striking two thin targets located near the beginning and the 
end of the fiducial decay region, producing 7r° and r] with known decay positions. From 
two-photon decays of tc° and rj, the reconstructed vertex position can be computed using 
the 7T° or the rj mass value (the rj mass value is taken from ||20||), and compared to the 
nominal target positions. Continuum production of 2n° events is also used, with the ad- 
vantage that the final state is very similar to the one of kaon decays. The two targets give 
energy scales consistent to better than 10~ 4 . The uncertainty on the overall energy scale 
is estimated from these comparisons to be ±3 x 10~ 4 . The corresponding uncertainty on 
the double ratio is ±2 x 10~ 4 . 

Non-linearities in the energy response are studied using K e3 decays, where the elec- 
tron energy measured in the calorimeter can be compared to the momentum measured 
in the spectrometer, and using data from the r] runs. Parameterising the deviations from 
linearity as AE/E — oc/E ± [3 • E, a is constrained to ±10 MeV and (3 is bound to be 
in the range ±2 x 10 -5 GeV -1 . Taking also into account larger deviations from linear- 
ity observed in the region E 1 < 6 GeV, the resulting uncertainty on the double ratio is 
±3.8 x 10" 4 . 

The uniformity of the calorimeter response over its surface is optimised using K e 3 
decays and checked using 7r° decays from the rj runs. Bias on the double ratio can arise 
from a dependence of the energy response on the photon impact radius r. Parameterising 
this effect as AE/E = 7 ■ r, 7 can be bound to be in the range ±10 -3 m _1 . In the 
region close to the beam tube, residual variations are smaller than 0.2%. The systematic 
uncertainty on the double ratio from these effects is ±1.6 x 10 -4 . 

Uncertainties in the correction of energy leakage from one cluster to another can lead 
to an apparent non-linearity and bias the double ratio. The correction used in the data 
is based on the transverse shower-profile measured during special runs in which single 
monochromatic electrons were sent to the calorimeter. The uncertainty in the shower 
profile is taken to be the difference between this measurement and the prediction of 
the GEANT Monte Carlo simulation. The resulting uncertainty on the double ratio is 
±1.1 x 10" 4 . 

The measurements of photon positions and the transverse size-scale of the calorime- 
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ter are adjusted and checked using K e3 decays, comparing the reconstructed cluster posi- 
tion with the electron track impact point extrapolated to the calorimeter. The associated 
uncertainty on the double ratio is ±1.6 x 10~ 4 . 

In the computation of the decay vertex position, the photon positions must be ex- 
trapolated to the longitudinal position of the maximum of the shower to account correctly 
for deviations of the photon directions from the projectivity of the calorimeter. Compar- 
ing data and Monte Carlo simulation in K e % decays, the uncertainty on this position is 
±2 cm. The resulting uncertainty on the double ratio is ±1.6 x 10 -4 . 

Finally, the effect of non-Gaussian tails in the energy response is minimised by 
the choice of the procedure used to adjust the overall energy scale. Residual effects are 
investigated by applying to the Monte Carlo samples a parameterisation of non-Gaussian 
tails in the energy response (arising mostly from photo-production of hadrons early in 
the electromagnetic shower) derived from K e s and rj data. No bias on the double ratio is 
observed within the Monte Carlo statistical error which is ±1.0 x 10~ 4 . 

Adding all the above uncertainties in quadrature, the total systematic error on the 
double ratio from the measurements of the photon energies and positions is found to be 
±5.3 x TO" 4 . 

For 7r + 7r~ decays, the vertex position is measured from the reconstructed tracks 
and is completely determined by the detector geometry. As a check, the reconstructed 
anticounter position can be measured in K$ — > 7r + 7r~ decays. The value obtained agrees 
with the nominal position to better than 1 cm Uncertainties in the geometry of the 
detector are at the level of 2 mm for the distance between the first two drift chambers, 
and 20 /xm/m for their relative transverse scale. This corresponds to possible deviations 
of « ±2 cm on the reconstructed AKS position. The corresponding uncertainty on the 
double ratio measurement is ±2.0 x 10~ 4 . The asymmetry in the Ks and Ki event losses, 
which could arise from the effect of non-Gaussian tails in the pt' 2 resolution, is smaller 
than 2.0 x 10 -4 . The overall uncertainty on the double ratio from the reconstruction of 
7r + 7r~ decays is therefore ±2.8 x 10~ 4 . 

6 Intensity effects 

6.1 Uncertainty on R due to accidental effects 

Most of the accidental activity in the detector is related to kaon decays in the 
high-intensity Ki beam. The overlap of extra particles with a good event may result 
in the loss of the event in the reconstruction or the selection Q, depending on the time 
and space separation of the activity in the detector. This effect is minimised by the 
simultaneous collection of data in the four channels and by the fact that Ks and weighted 
Kl decays illuminate the detector in a similar way. The possible residual effect on the 
double ratio can be separated into two components: intensity variations between the two 
beams coupled to different, intensity-dependent, event losses in the n + ir~ and 7r°7r° modes 
(intensity difference), and a residual difference in the illumination between Ks and Kl 
decays coupled to a variation of the event loss with the impact points of the K° decay 
products (illumination difference). 



6) The difference was 2 cm in the 1998-99 data. 

7 ) The fosses in the trigger are already accounted for in the measurement of the trigger efficiency. 



11 



6.1.1 Intensity difference 

If the losses depend, as expected, linearly on the Kl beam intensity, the intensity 
difference effect is given by: 

AR = AX x AI/I (5) 

where AA is the difference between the mean losses in tt + tt~ and 7r°7r°, and AI/I is the 
difference in the mean Kl beam intensity as seen by Kl and K$ events. 

The accidental rate can be measured directly from the activity in the detector 
within the readout time-window before each event. Comparing the rate of out-of-time 
LKr clusters and out-of-time tracks in good K$, Kl — > 7r + 7r~ decays, AI/I is found to be 
respectively (+0.4 ± 0.4)% and (+0.6 ± 0.3)%, where the quoted uncertainties are only 
statistical. For this measurement Kg and Kl are identified by the decay vertex position 
to avoid the correlation between mistagging probability and beam intensity. The bias on 
AI/I from the fact that this measurement uses good reconstructed decays is a negligible 
second-order effect. The measurement of accidental activity is illustrated in Figure [| 
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Figure 2: Probability of accidental activity in the LKr in the ~ 150 ns readout window, 
as a function of the time during the spill, for Kg and Kl decays to 7r + 7r~. 

The Kl beam intensity for each event can also be estimated using the information 
from the Kl beam monitor, integrating the intensity over a 200 ns time window. The 
difference between the average intensities as seen by K$ and Kl decays is found to be 
(—0.08 ± 0.04)%, where the quoted uncertainty is only statistical. The systematic uncer- 
tainty on this measurement is estimated to be less than 1%. Another method to study 
the correlation of the two beams, which does not rely on the use of good events, consists 
of a direct computation of the correlation between the Kl and K$ beam monitor counts, 
using events taken uniformly in time. To avoid statistical fluctuations in the rate mea- 
surements, this is done using a 15 /is integration time of the beam monitors. All known 
Kl/Ks variations take place on a much longer time scale. This method confirms that 
AI/I is consistent with within 1%. The final estimate for AI/I is (0 ± 1)%. 

The beam-induced event losses can be evaluated by overlaying data with beam mon- 
itor (BM) triggers taken in proportion to the beam intensity, which reproduce the ambient 
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activity as seen by the detectors. Detector information from the original events and the 
BM triggers are superimposed at the raw-data level and the overlayed events undergo 
the standard reconstruction and analysis procedure. Losses and gains from migration of 
events around the cuts are taken into account. When estimating the event losses related 
to accidental activity, the losses related to doubling the noise effect inherent in this pro- 
cedure must be removed. The net event losses induced by the accidental activity are at 
the level of 1 to 2%. The dominant source of event loss is found to be the appearance of 
the drift chamber overflow condition in the overlayed event. This loss is higher for tt + tt~ 
decays than for 7r°7r° because of the presence of at least two drift chamber hits per plane 
for the original tt + 7t~ decay. From this procedure, AA is found to be 1.0%. The overlay 
procedure can also be applied starting from a Monte Carlo sample of original 7r + n~ and 
7t°7t° decays. In this case, AA is found to be 0.65%. The statistical uncertainties on AA 
from the overlay samples are less than 0.1%. The difference between the two estimates is 
attributed to the lower hit-multiplicity in the Monte Carlo samples compared to the data. 

AA can also be estimated directly from the data, by comparing the ratio of 7r + 7r~ 
and events obtained in the normal Kl + Ks beam runs and in pure K$ runs, in which 
the accidental activity is one order of magnitude smaller. This leads to AA = (0.9±0.6)%. 
Similarly, AA can be checked by dividing the data into bins of Ki beam intensity and 
looking at the variations in the ratio of 7r + 7r~ and 7r°7r° events. This leads to an estimate 
in agreement with the values above. In conclusion, from this study the final estimate of 
AA is (1.0 ± 0.5)%. This value is typically 30% lower than the estimate derived for the 
1998-99 data-taking period, as expected from the lower beam intensity in 2001. 

Finally, the linearity of the losses with intensity can be checked by looking at the 
losses as a function of the beam intensity given by the BM triggers. Figure |3] shows the 
accidental event losses in the ti + it~ and 7r°7r° modes as a function of the beam intensity 
from the overlay procedure applied to Monte Carlo events. 

Taking into account all the above results, the estimate of the uncertainty on the 
double ratio related to differences in intensity-dependent losses is ±1.1 x 10~ 4 

6.1.2 Illumination difference 

The illumination-difference effect has been estimated from the overlay samples, com- 
puting separately the losses for Ks and Kl events. This computation has been performed 
using both data and Monte Carlo original events. In the first case, the value obtained is 
(+0.9±3.5) x 10~ 4 , in the second it is (±1.4±2.8) xl0~ 4 , where the quoted uncertainties 
are the statistical errors from the overlay samples. As expected, there is no evidence of a 
significant effect, and we use as uncertainty on the double ratio ±3.0 x 10~ 4 . 

6.1.3 Overall uncertainty on R and cross-check 

Combining the two above uncertainties in quadrature, the total uncertainty on R 
from accidental effects is ±3.1 x 10~ 4 . This uncertainty is dominated by the statistical 
error of the overlay procedure. 

The overlay method can also be used to estimate the combination of the two effects. 
For this, Kg events are overlayed only with BM triggers from the Ks monitor and K^ 
events only with BM triggers from the Kl monitor. This method relies on the accuracy of 
the BM triggers to estimate correctly the Kl-Ks intensity difference. From this method, 
the overall accidental effect on the double ratio is found to be (4.7 ± 4.9) x 10~ 4 , where 

8 ) For the 1998-99 sample, this uncertainty was estimated to be ±3.0 x 10~ 4 . 
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Figure 3: Beam-induced net event losses as a function of instantaneous intensity in arbi- 
trary units. Lines are drawn to guide the eyes. 

the quoted error is the statistical uncertainty from the overlay sample. The systematic 
uncertainty from the accuracy of the BM triggers is expected to be < 2xl0~ 4 . From this 
cross-check, there is no evidence of any unexpected effect on the double ratio. 

6.I.4 In-time activity from the K$ target 

The techniques above do not take into account any additional detector activity in 
the Ks beam generated by the same proton which produced the K$ event. Studies of this 
background, mostly searching in the LKr calorimeter for additional clusters in 27r° events 
from pure K s beam runs, allow an upper bound on the effect on R of 1 x 1CT 4 to be set. 

6.2 Origin of Aa LS 

The accidental tagging probability ais depends only on the proton beam intensity 
seen by the Tagger and, consequently, it is to first order the same for 7r°7r° and 7t + tt~. 
However, because the event selection is more sensitive to accidentals for 7r + 7r~ events, we 
expect a difference of the measured oils f° r the 7r°7r° and 7r + 7r~ samples due to beam- 
intensity variations with time. A quantitative understanding of the effect of accidental 
activity on selected events can be reached by studying the BM trigger overlays. 

A value Aais — (+3.5 ±0.4) x 10 -4 is expected from the overlay of data, where the 
error is only statistical. Using loss and gain probabilities from the overlay Monte Carlo 
samples instead of data we find Aais = (+2.0 ± 0.4) x 10 -4 . 

In Figure |] the variation of Aa^s within the spill is shown and compared with the 
overlay computation. Most of the difference in accidental tagging between the 7r°7r° and 
7r + 7r~ modes comes from the beginning of the spill where the instantaneous intensity is 
higher and the beam intensity variations more significant. 

Another source of event losses is the inefficiency of the 7r + 7r~ trigger. The intensity- 
dependent part is studied separately and its effect on Ao>ls is estimated to be (0.4 ± 
0.2) x 10- 4 . 
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Figure 4: Measured compared to predicted values of Aais as a function of the time during 
the spill 

The observed Aoils value of (+3.4 ± 1.4) x 10~ 4 is therefore well reproduced both 
qualitatively and quantitatively. 

7 Result 

The effect on the result of the corrections described above, and the various sources 
of systematic uncertainties are summarised in Table [l]. 



Table 1: Corrections and systematic uncertainties on the double ratio R (2001 data) 





in 10" 4 


7r + 7r trigger inefficiency 


+5.2 


± 3.6 


(stat) 


AKS inefficiency 


+1.2 


± 0.3 




Reconstruction °^ n + 

Of TX^TX 




± 5.3 
± 2.8 




Background ^° n 7 

tO 7T + 7T 


-5.6 
+ 14.2 


± 2.0 
± 3.0 




Beam scattering 


-8.8 


± 2.0 




Accidental tagging 


+6.9 


± 2.8 


(stat) 


Tagging inefficiency 




± 3.0 




. statistical 
Acceptance 

systematic 


+21.9 


± 3.5 
± 4.0 


(stat) 


. ., , , ; . . t intensity difference 
Accidental activity .„ . 

illumination difference 




± 1.1 
± 3.0 


(stat) 


Ks in time activity 




± 1.0 




Total 


+35.0 


± 11.0 





Figure |5] shows the measured double ratio after corrections as a function of the kaon 



energy. 
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Figure 5: Measured double ratio R in kaon energy bins. 

The final result for the double ratio from the 2001 data set is R = 0.99181 ± 
0.00147 ± 0.00110, where the first error is the statistical error from the 27r samples, and 
the second is systematic. Out of this systematic uncertainty, ±0.00065 is due to the finite 
statistics of the control samples used to study the systematic effects. 

Many cross-checks of the stability of the result have been performed, by varying 
some of the selection cuts and by searching for a dependence of the result on several 
variables, such as the beam intensity, the time during the spill, and the data-taking period. 
No significant variation in the result is observed. An analysis adopting a different scheme 
for data compaction, filtering, selection and correction procedures was also performed in 
addition to the one presented here. Its result fully confirms the above measurement. 

The corresponding value of the direct CP- violation parameter Re(e'/e) from For- 
mula |l] is: 

Re(e'/e) = (13.7 ± 2.5 ± 1.1 ± 1.5) x 10" 4 , 

where the first uncertainty is the pure statistical error from the 2ir samples, the second is 
the systematic error coming from the statistics of the control samples, and the third is the 
contribution of the other systematic uncertainties. Combining the errors in quadrature, 
the result is: 

Re(e'/e) = (13.7 ±3.1) x 10~ 4 . 

This result is in good agreement with the published value from the 1997-98-99 data: 
Re(e'/e) = (15.3 ± 2.6) x 10" 4 . 

The comparison of the present and earlier results is particularly significant since they 
were obtained from data taken at different average beam intensities. The correlated sys- 
tematic uncertainty between the two results is estimated to be ±1.4x 10 -4 . Taking this cor- 
relation into account, the combined, final result on Re(e'/e) from the NA48 experiment is: 

Re(e'/e) = (14.7 ± 1.4 ± 0.9 ± 1.5) x 10" 4 , 

or, with combined errors: 

Re(e'/e) = (14.7 ±2.2) x 10~ 4 . 
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